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ABSTRACT 



Context. The recently-discovered lack of close binaries, among extreme horizontal branch (EHB) stars in Galactic globular clusters, 
has thus far constituted a major puzzle, in view of the fact that blue subdwarf stars - the field counterparts of cluster EHB stars ~ are 
well-known to present a high binary fraction. 

Aims. In this Letter, we provide new results that confirm the lack of close EHB binaries in globular clusters, and present a first scenario 
to explain the difi'erence between field and cluster EHB stars. 

Mettiods. First, in order to confirm that the lack of EHB binaries in globular clusters is a statistically robust result, we undertook a 
new analysis of 145 horizontal branch stars in NGC6752, out of which forty-one belong to the EHB. To search for radial- velocity 
variations as a function of time, we repeated high-resolution (R = IS 500) spectroscopy of all stars, four times during a single night 
of observations. 

Results. We detected a single, hot (25 000 K), radial-velocity variable star as a close-binary candidate. From these results, we estimate 
an upper-limit for the close (period P < 5 day) binary fraction / among NGC6752 EHB stars of 16% (95% confidence level), with 
the most probable value being / = 4%. Thus our results clearly confirm the lack of close binaries among the hot HB stars in this 
cluster. 

Conclusions. We suggest that the confirmed discrepancy between the binary fractions for field and cluster EHB stars is the con- 
sequence of an /-age relation, with close binaries being more likely in the case of younger systems. We analyze theoretical and 
observational results available in the literature, which support this scenario. If so, an age difference between the EHB progenitors 
in the field and in clusters, the former being younger (on average) by up to several Gyr, would naturally account for the startling 
differences in binary fraction between the two populations. 
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1. Introduction 



The presence of a large population of binaries among field 
B-type subdwarf (sdB) stars, also referred to as Extreme 
Horizon tal Branch (EHB) sta rs, is w ell-establis hed in the lit- 
eratur e dFerguson et al. 198 ^ lAllard et al. 19 91 lUHa & Theilll 



1998'; 'Aznar Cuadrado & Jefferv ToOl; Maxtedetal.' '2001' 
Williams et al. 2001; Reed & Stiening 2004; Napiwotzki et al. 
20041) . The measurement of binary fraction varies from one sur- 



vey to another, probably because of the different survey selection 
effects. It is generally agreed however, that the binary fraction 
must be large. Moreover, extensive analysis o f orbital param- 
eters (Moraneta l. 1999; Safferetal. 1998; Hebe r et al.| [ 20021; 
iMorales-Rueda et al. 2004, 2006) reveal that close binaries (P < 
10 d) comprise a major fraction of the field EHB stars, with 
short-period systems — including either a degenerate or a low- 
mass main sequence (MS) companion — constituting about half 
of the entire field sdB population. Indeed, theoretical results in- 
dicate that these stars can be naturally explained wit hin the con- 
text of binary star evolution jHan et alj|2002ll2003b . Moreover, 



* Based on observations with the ESO Very Large Telescope at 
Paranal Observatory, Chile (proposal ID 075.D-0492(A)) 



EHB stars of binary origin may also account for th e UV flux 
excess ("UV upturn") observed in elliptical galaxies dHan et al.l 
120071). although it shoul d be noted that single-star scenarios have 
also been suggested (see lCatelanll2007L for a recent review). One 
way or another, the link between field sdB stars and binary sys- 
tems is very well-established, both on observational and theoret- 
ical grounds, and close, short-period systems are the most fre- 
quently found amongst them. 

It therefore came as a great surprise that the first radial- 
velocity (RV) surveys among EHB stars in globular clusters 
(GCs) revealed a remark able lack of close binary systems (see 
Im oni Bidin et al.ll2007al for a review). 

In this Letter we present, in light of new observational re- 
sults, a first scenario to account for this discrepancy. 



2. Observations and measurements 

High-resolution spectra {R - 18 500) centered on H^ (4700- 
4970 A) were collected with 3600s exposures on June 29, 2005 
at VLT-UT2, with the FLAMES -GIRAFFE optical spectrograph 
in MEDUSA configur ation (setup HR7A ). 169 HB stars were 
selected from the iMomany et al.l (l2002h photometry, spanning 
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Fig. 1. Results of our RV variability search. The maximum RV variation (RVmax - RVmin) observed for each star is plotted at its 
effective temperature. For the sake of clarity, we plot only the 3cr error bars, which indicate the statistical significance of the observed 
variations. The data for the only star with a clearly-detected RV variation, is plotted using a larger, solid circle symbol, and thicker 
lines to indicate its errorbar. 



Table 1. Time elapsed (in hours) between the start of each ex- 
posure and the first one, for each fiber configuration (Medusa 1 
and Medusa 2). 



frame 


Medusa! Medusa2 




hours 


1 


0.00 0.00 


2 


1.03 4.73 


3 


4.75 5.77 


4 


5.78 9.13 



the whole temperature range along the HB, out of which 54 
were EHB stars (i.e., with Teff > 20000 K). When perform- 
ing the target selection, we ensured that each individual candi- 
date was sufficiently isolated that no other bright star fell within 
the fiber, extending observations as far as possible toward the 
cluster center. Selected targets were divided into two fiber con- 
figurations. In Table [T] we indicate the time at the start of each 
exposure, in hours from the first spectrum. To confirm that our 
targets were true cluster HB stars, we verified that absolute ra- 
dial velocities are consistent with cluster membership, and we 
inspected all spectra by eye. We excluded six, cool-star targets 
(Teff ~ 6 500 K) from our analysis, because their cluster mem- 
bership was doubtful. 

All data reduction steps were completed using ESO's 
GIRAFFE pipeline. Wavelength calibration was performed us- 
ing lamp fiber data acquired simultaneously during the observa- 
tions: according to the FLAMES manual the calibration error is 
about 0.15 km s"'. 

RV 's were measured usin g the cross-correlation (CC) tech- 
nique (' Tonrv&DavisllI979l) . cross-correlating each spectrum 
with synthetic templa tes of approp riate temperature and grav- 
ity from the library of iMunari et al. 1 12005). During CC we used 
only H/3 with full wings, neglecting all weaker lines that were 



hardly visible due to noise, and highly variable from star to star 
Some tests carried out prior to the actual measurements revealed 
that, with this limitation, the adopted metallicity of the template 
had a negligible impact on the results. Errors in the measured 
RV variations were calculated as the quadratic sum of the errors 
for the individual measurements. 

More details about the observations, data reduction and mea- 
surements will be provided in a forthcoming paper (Moni Bidin 
et al., in preparation). 

3. Results 

We were unable to acquire observations for twelve stars, due 
to for example, light contamination by nearby fibers, too low 
signal-to-noise ratio, or the star being badly-centered in the fiber 
Moreover, for six of the hotter stars one or more spectra were of 
too low S/N to provide reliable measurements; the observations 
being thus incomplete for these stars, we excluded them from 
subsequent analysis. Therefore, we successfully performed a RV 
analysis for 145 stars, out of which 41 fall along the EHB. 

Our analysis revealed that the errors given by CC alone 
were underestimated. We evaluated this effect by analyzing sta- 
tistically the distribution of for the cooler stars (Tetf < 
10 000 K). Among this subsample, we do not expect a population 
of very close binaries, and observed variations can be ascribed 
only to random errors. We found that the latter were underesti- 
mated by a factor of 1 .44. We expected this factor to increase for 
decreasing S/N, but repeating the analysis in different tempera- 
ture bins we found only small variations. Moreover, evaluating 
the correction factor for hotter stars would invalidate our results, 
because RV variations due to binarity, if not properly identified, 
would spuriously enlarge the estimated errors and thus incor- 
rectly lead to the lack of detected binaries. Hence we applied 
the correction by a factor 1.44 to the whole sample, as a good 
approximation valid for all temperatures. 
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Fig. 2. Probability of detecting a close binary systems in our sur- 
vey, as a function of its period. As probability varies with both 
times of observations and errors in the measured RV variations, 
we plot both the average over the probabilities for each individ- 
ual star (solid line) and the one calculated assuming a fixed error 
of 10 km s"' for all stars (dotted line). 



Our results are summarized in Figure [T] For each star, we 
plot the maximum RV variation and the 3cr interval. Effective 
te mperatures were e s timated using the color- Teif relation derived 
in iJMoni Bidin et all (l2006bh . Because each datum is the maxi- 
mum of six measurements, the probability that each datum ex- 
ceeds 3cr is 1.6%, and we expect 0.66 false detections out of 41 
EHB targets. We adopted this threshold to select objects display- 
ing suspect variations to be further analyzed, but we found no 
need to discuss border cases because the only candidate found 
shows variation at the lOcr level (16.1 + 1.6 km s"'). This star 
was akeady highlighted by Moni Bidin et al. (2007b) for its red 
color and evident Mg lb triplet, signatures of a cool MS com- 
panion. We are currently analyzing data for this star, which turns 
out to be the first sdB+MS close system ever found in a globu- 
lar cluster In summary, we come out with only one good close 
binary candidate in our whole sample of 145 HB (41 EHB) stars. 

In Figure |2] we plot the estimated detection probability in 
our survey a s a func tion of the period. This was calculated as 
m oni BidinetaD (^006b), assuming a mass of 0.5 Mq for 
both components and circular orbits, and integrating over the 
phase and the orbital inclination angle. We derived the detec- 
tion probability for each individual star (given sampling times 
in Table [T] and the aforementioned 3cr values), and then aver- 
aged it to obtain the final curve. We also checked that other rea- 
sonable ways of tackling the problem (for example, assuming a 
fixed 3cr = 10 km s ' for all stars in the sample) did not change 
our results by more than 1-3%. Clearly, the detection probability 
rapidly falls down for increasing periods, as a consequence of 
the observations being restricted to a single night; accordingly, 
we will limit our analysis to systems with P < 5 d. It must be 
emphasized that such a limitation is not particularly severe, since 



many surveys have revealed that the bulk of the field binary sdB 
population has periods of approximately f ^ 1 d, systems with 
P > 5 d constitutin g only the tail end of the period distribution 
(see, for example, Morales-Rueda et aLll2003b . 

Using the calculated detection probability, we estimated the 
close binary fraction /p<5 d among EHB stars as in lMaxted et al.l 
(1200 1), assuming a flat distribution of periods, and the successful 
detection of one binary out of 41 targets. We cal culated /f<5d 
assum ing a Gaussian distribution i n log P, as in iMaxted et al.l 
(1200 1') and'N apiwotzki et all (|2004 . which better represents the 
period distribution for field sdBs — but the results do not differ 
noticeably. The derived probability distribution peaks at / = 4%. 
It is markedly non-Gaussian, but it falls below 5% for / > 16%. 
We conclude that the best estimate (i.e., the most likely value) 
for fp<5d is 4%, and that /p<5d <16% at the 95% confidence 
level. Our derived upper limit is in perfect agreement w ith the 
one that was found for //><iod bv iMoni Bidin et akl (l2006b). We 
note that these results a re at variance with th e preliminary ones 
previously obtained by iPeterson et al.l ( |2002|) . who claimed the 
detection of many close binaries among the EHB stars in the 
same cluster The reader is referred to Moni Bidin et al.. (20063) 
for a discussion of this discrepancy. 

4. Is there a close binary fraction-age relation? 

Usin g an independent dataset, a sample more than twice as large 
as in iMoni Bidin et al. I (I2006bh . and a resolution in RV varia- 
tions that is higher by almost a factor of two, for the first time 
we were able to find a good binary candidate among the EHB 
stars in NGC 6752. That notwithstanding, our results confirm 
that the corresponding (close) binary fraction / is very small 
in NGC 6752, with a most likely value of / = 4% and an upper 
limit of f - 16% at the 95% confidence level. 

There are hints that a small / is not a peculiarity of 
NGC 67 52, but could also be a c haracteristic of other globular 
clusters 

dM 

oni Bidin et al ] |2006ah . This is in sharp contrast with 
the situation for field sdB stars, where close binaries are at least 
a factor of ten more frequent, comprising up to 70% of the en- 
tire sdB population (see §4 in Catelan 2007, for a recent review). 
There is however no knowledge about cluster EHB stars in long- 
period binaries, or with a close low-mass companion, because 
no survey has investigated their role yet. These kind of systems 
are known to exist among field sdBs, but are just a minor popu- 
lation, and their presence in GCs would not alleviate the striking 
contrast with field results. 

What is the origin of this startling difference between field 
and cluster EHB stars? We believe that it may not be completely 
unexpected: there should be a relation between the close bi- 
nary fraction and the mean age of a sdB population, as a conse- 
quence of the different efficiency of binary channels responsible 
for EHB star formation. 

Theoretical arguments strongly suggest that sdB stars in 
close binary systems should have u ndergone a t leas t one com- 
mon envelope (CE) phase. Although lHan et akl (I200Z) found an 
upper limit for the initial mass of the sdB progenitor in this sce- 
nario, they also pointed out that within the permitted values a 
higher initial mass favors the CE channel, and leads to sdB bina- 
ries with shorter periods. In fact, a higher mass implies a more 
tightly-bound envelope, which requires a great er amount of (or- 
bital) energy to be released. On the other hand. lHan et al.l(l2002l) 
explored the stable Roche Lobe Overflow (RLOF) scenario, and 
found that a higher initial progenitor mass makes it harder for 
the RLOF to be stable (see their Table 3), because the mini- 
mum mass of the companion increases with increasing progeni- 
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tor mass. For higher values, fewer MS and white dwarf (WD) 
secondaries are sufficiently massive for the activation of this 
channel (sdB's with neutron star companions are indeed very 
rare). The progeny of systems that underwent stable RLOF, that 
is wide binaries with very long periods, is therefore generated 
primarily by progenitors of lower initial mass. 

In light of these results, a relation between / and the mean 
initial mass of the sdB progenitors could be naturally expected, 
hence implying a relationship between / and (mean) age. More 
specifically, one may naturally expect that field sdB's formed 
from progenitors with a wide spectrum of initial masses (up to 
about 2 Mq), whereas in an old population (such as in globu- 
lar clusters) only the progeny of less massive stars are currently 
found on the EHB, those of more massive ones having long 
evolved away from the He-burning phase. The stable RLOF is 
an efficient channel for sdB formation in the old case, while the 
CE one is not — and the CE itself would be released at earlier 
stages, before the orbits shrink substantially. Therefore, in old 
populations we should expect to find predominantly wide bina- 
ries, or/and single EHB stars formed through other channels (see 
[Catelan 2007, for a recent disc ussion). In fact, WD mergers, the 
third binary channel studied bv lHan et al.l ( l2002h . can form (sin- 
gle) sdB stars, and is expected to be particularly efficient for old 
populations (see below). Moreover, it may be expected to play an 
important role within the dense environment of a globular clus- 
ter, w here stellar encounters can harden close binaries (Heggie 
Il975h and thus enhance channel efficiency, as proposed earlier 
by Bailvn & Ibenl(ll989l) . However, dyn amical effects are no t re- 
quired for the formation of EHB stars dWhitnev et al.|[T998l) . as 
indicated by t he lack of ra dial gradients among the EHB stars 
in NGC2808 jBedin et al.| [2000) and w Centauri (NGC5139; 
lD'Cruzetal.ll2000l) . 

There are other results in the recent literatu re that pro- 
vide su pport to our proposed scenario. In fact, Napi wotzki et alj 
(|2004|) already invoked a possible /-age or /-metallicity rela- 
tion to explain their results. Among field sdB stars, they found a 
close binary fraction lower than in previous surveys, and noted 
that their sample was on average fainter, possibly including more 
thick disk or halo (hence older/metal poorer) stars. This sug- 
gested difference between the data used has never been fully in- 
vestigated, and a further check of the populations of these stars 
would be invaluable. Even stroriger sup port is provided by the 
recent simulations by Han et all ( 1200 7). modeling the UV up- 
turn of elliptical galaxies from a binary sdB population. Their 
Figure 7, though aimed at reproducing the spectral energy distri- 
bution of elliptical galaxies, shows how the relative contribution 
of sdB's formed through different channels evolves with time. 
One finds that the contribution to the UV flux from RLOF sdB's 
is essentially constant with time, whereas the one from sdB's that 
experienced the CE phase first becomes important for a popula- 
tion of age 1 .5 Gyr and then slowly fades, being more and more 
marginal for increasing ages. On the other hand, for populations 
older than 5 Gyr these authors find that sdB's formed through 
the merger channel dominate — ^but these are not binaries any 
longer. We suggest that these theoretical results implicitly hint 
at a solution to the heretofore puzzling lack of close EHB bina- 
ries in globular clusters, on the one hand, and their large numbers 
among field stars, on the other 

As an alternative to the binary scenario, it is worth noting 
that many single-star evolutionary channels have been invoked 
to explain EHB star formation in GCs, including interactions 
with a close planet (Soker 1998, see also Silvotti et al. 2007), 
He mixing drive n by internal rotation (Sweigart & Mengel 1979; 
ISweigarll Il997h or by stellar encounters dSuda et al.. ,2007i) . 



dredg e-up induce d by H-shell instabilities ( yon Rudloff et al.l 
Il988l but see alsol Denissenkov & VandenBergll2063h . close en- 
count ers with a central, intermediate-mass black hole dMiocchil 
l2007h. and a sub-population of stars with high helium abundance 
(e.g.. lD'Antona et alj|2005l) . 

We note that any EHB formation mechanism involving sin- 
gle progenitors should be particularly inefficient among field 
stars (but possibly not so among GCs), since the EHB progeni- 
tors are on average younger, and thus more massive, in the field 
than in GCs. As a consequence, a much higher amount of RGB 
mass loss is needed to successfully produce a EHB star in the 
case of a field progenitor, whereas the ancient GC red giants 
have a much smaller envelope to be removed before they be- 
come EHB stars. This notwithstanding, a small component of 
single-star progeny may still be ne eded among field st ars to fully 
explain the available observations jLisker et al.ll2005l) . 
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